Prediction of Composite Laminate Strength Properties Using 

a Refined Zigzag Plate Element 


A. Barut 1 , E. Madenci 2 

University of Arizona, Tucson, AZ, 85721, USA 
and 

A. Tessler 3 , 

NASA Langley Research Center, Hampton, VA 23681-2199, USA 


This study presents an approach that uses the refined zigzag element, RZE (2,2) in 
conjunction with progressive failure criteria to predict the ultimate strength of composite 
laminates based on only ply-level strength properties. The methodology involves four major 
steps: (1) Determination of accurate stress and strain fields under complex loading 
conditions using RZE 2 21 -based finite element analysis, (2) Determination of failure locations 
and failure modes using the commonly accepted Hashin’s failure criteria, (3) Recursive 
degradation of the material stiffness, and (4) Non-linear incremental finite element analysis 
to obtain stress redistribution until global failure. The validity of this approach is 
established by considering the published test data and predictions for (1) strength of 
laminates under various off-axis loading, (2) strength of laminates with a hole under 
compression, and (3) strength of laminates with a hole under tension. 


I. Introduction 

Laminated fiber-reinforced polymer (FRP) composites exhibit directional stiffness and strength properties, and 
offer different fabrication architecture; thus, enabling advanced design concepts, structural tailoring, multi- 
functional features, and performance enhancements. Their actual applications, however, are often lagging behind 
the material development and fabrication because characterization of their mechanical and strength properties 
requires large amounts of data to achieve an adequate confidence level. The ability to model and simulate 
composite materials for strength prediction based on ply level material property data can reduce the time and cost 
for material testing and characterization; thus, provide the opportunity to use advanced material systems on critical 
military and commercial equipment in a timely manner. 

Strength of FRP composite laminates depends on fiber directions, type of fibers and resins, number of layers, 
stacking sequence, type of loading, environmental conditions, and random variation of material and strength 
properties due to manufacturing. Their failure involves a progressive series of events with discrete failure modes 
such as matrix cracking, fiber-matrix shear, fiber breakage, and delamination. The presence of such failure modes 
results in stiffness degradation. This leads to stress redistribution in the layers and constituents. Therefore, to enable 
realistic simulation of composite -laminate strength, accurate predictions of stress redistribution is paramount. 

The stress state in a composite laminate is dependent on the loading conditions, fiber- and resin-rich layer 
thickness, face-sheet lay-up, core thickness and process-dependent properties. Understanding the behavior of such 
structures can result in weight and cost savings by designing against unnecessary conservatism. Finite element 
Analysis (FEA) is widely utilized to assess the strength of such structures. To model a vast majority of the potential 
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failure modes, however, an extremely detailed finite element discretization is often required. Using traditional 
elements, through-the -thickness finite element discretizations are often impractical. This is because an extremely 
high-density mesh is required to maintain a proper aspect ratio between the elements in the fiber- and resin-rich 
layers. Typically, composite laminates have a low ratio of thickness to in-plane characteristic dimension; thus, they 
can be modeled using 2-D plate elements that are based on an assumed set of kinematical relations. However, they 
require shear correction in order improve transverse shear deformation. The constitutive relations are in terms of 
average stress resultant quantities rather than point stresses and strains. Strength prediction requires an accurate 3-D 
stress field especially if edge delamination is the source of progressive failure. 

Material characterization of laminated composite structures is well established, and failure prediction models 
for damage modes associated with FRP composite laminates are relatively successful. Ochoa and Reddy [6] present 
an excellent discussion of progressive failure analyses. In laminated composite structures, fiber-dominated failure 
modes (e.g., axial tension) are better-understood and analyzed than matrix-dominated failure modes (e.g., axial 
compression, transverse tension). To predict ultimate strength of a laminate, there exist two commonly accepted 
methods: damage mechanics and progressive ply failure. The damage-mechanics approach employs physically 
based equations for damage initiation and evolution while taking into account the material’s microstructure [1-3]. 
However, it requires extensive material characterization for damage parameters, and in most cases such 
measurements are not feasible. Progressive ply failure combines failure criteria for damage mode identification [4- 
5] and degradation [6-9] of the material stiffness. The value of degradation factor is not-physically based; it is 
assigned a small enough value in order to ensure the convergence of finite element analysis. Methods to predict 
failure initiation and to perform material degradation remain active areas of research. Material degradation can be 
performed using a ply-discounting approach or an internal state variable approach based on continuum damage 
mechanics. Most composite failure analysis methods embedded within a finite element analysis tool perform a 
point-stress analysis, evaluate failure criteria, possibly degrade material properties, and then continue to the next 
solution increment. 

There exist many progressive failure models in the literature to predict the strength of FRP composite materials. 
These previous models utilize ply level strength parameters; however, they lack reliability in their strength 
predictions under longitudinal and transverse tension and compression, longitudinal and transverse open-hole 
tension and open-hole compression. The lack of reliability is associated with inaccurate computation of stress 
redistribution within the conventional finite elements used to model composite laminates. The existing conventional 
elements based on the first order shear deformation theory (FSDT) require different shear correction factors in the 
presence of a degraded ply or plies through the thickness, and fail to provide accurate and complete stress field near 
the free edges. 

To evaluate the strength of a composite structure with a high degree of fidelity, it is first essential to predict 
accurate stress and strain fields, and then utilize a computationally efficient and accurate progressive damage model. 
For these reasons, an accurate and robust finite element analysis is required that can account for the discrete nature 
of both the fiber- and resin- rich layers of the individual plies, as well as the variation of stiffness and strength 
properties of the core. Recently, a viable theoretical framework called the Refined Zigzag Theory (RZT) was 
developed by Tessler et al. [10-13], The theory uses a fixed number of kinematic variables (seven for plate analysis) 
and is well suited for the analysis of a wide range of composite and sandwich laminates without the use of shear 
correction factors. Efficient RZT -based finite element formulations for the analysis of beam and plate structures 
were also developed [14-16], Focusing specifically on higher-order deformation effects commonly associated with 
thick-section laminated composite and sandwich panels, Barut et al. [17] explored a higher-order kinematic 
representation using RZT as its basis. In their higher -order RZT formulation, Barut et al. used the same zigzag 
functions as those of RZT; however, higher-order in-plane deformation modes were added using a set of piecewise - 
quadratic zigzag contributions and, moreover, both linear and quadratic polynomial expansions were added for the 
transverse displacement component. Hence, this higher-order RZT is defined by the notation RZT 12 21 , in which the 
first and second superscripts denote the order of expansions used for the in-plane and transverse displacement 
components, respectively. 

The focus of this study is to present a validated approach based on the refined zigzag element, RZE i2 21 by 
Barut et al. [17] used in conjunction with the progressive failure criteria of Hashin [4] to predict the ultimate 
strength of composite laminates utilizing the ply-level strength properties only. The validity of this approach is 
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established by considering the published test data and predictions for (1) strength of laminates under various off-axis 
loading [18, 19], (2) strength of laminates with a hole under compression [20], and (3) strength of laminate with a 
hole under tension,[8, 21-22], 


II. Approach 

This approach involves four major steps: (1) determination of accurate stress and strain fields under complex 
loading conditions using RZE 1221 based finite element analysis, (2) determination of failure site(s) and modes using 
commonly accepted failure criteria of Hashin, (3) recursive degradation of the material stiffness, and (4) non-linear 
incremental finite element analysis to obtain stress redistribution until global failure. 

Refined Zigzag Element, RZE 12 21 [17] 

This element offers substantial analytic and computational advantages for the analysis of homogeneous, 
laminated composite, and sandwich laminates. Starting from a multi-scale kinematic description represented by a 
superposition of the coarse and fine kinematics, it employs a parabolic transverse displacement assumption and 
higher-order RZT-based in-plane displacements to derive a higher-order refine zigzag plate theory for laminated 
composite and sandwich plates. The displacement assumptions involve eleven kinematic variables and are given as 


Mi (t) (x, z) = k(x) + zO l (x) + ^ (z) y/ x , (x) + y <fc k) (z) y/ n (x) 

h 


«2 (x. z) = V(x) + z9 2 (x) + (j)\ (z) y/ 2l (x) + (z) y 22 (x) 


h 


M 3 (x, z) = w(x) + - Wj (x) + 
h 




w 2 (x) 


( 1 ) 


where x = (x 1 ,x 2 ) are the in -plane coordinates, z e (—h, h) is the thickness coordinate, with 2 h denoting the total 
plate thickness; <p^\z) are the piecewise linear zigzag functions within the k' b layer that were originally derived in 
[10]. The kinematic variables m(x) and v(x) are the through-the -thickness uniform components of the in-plane 
displacements in the x- and y-dircctions, respectively. The kinematic variable w(x) is an average transverse 
displacement defined in [17]. The additional transverse-displacement variables, w i ( x ) and w 2 (x) , represent the 
symmetric and anti -symmetric thickness-stretch modes. The kinematic variables 9 t (x) and 0 2 (x) denote the 
bending rotations about the negative x and positive y axes, respectively. The zigzag amplitudes defined by the 
variables i// al (x) and y al (x) ( a = 1, 2 ) permit both non-symmetric, y/ al (x) , and symmetric, y/ a , (x) , in-plane 
zigzag deformation modes. The contributions associated with i// l2 (x) , i// 22 (x ) , w,(x), and w 2 (x) are the 
additional, higher-order terms that do not appear in the RZT model [10-13] . In addition, an average transverse 
normal stress, a , is independently assumed as a cubic function through the laminate thickness as 


CT.,(x, z) = cr „(x) + cr.(x) 


3 X 

z 

3h 3 , 


( 2 ) 


where a n (x) and a 2 (x) are functions of the kinematic variables defining the displacement field in Eq. (1) (refer to 
[17] for the details of derivation of these functions.) 

The zigzag kinematic assumptions are represented by piecewise-linear C°-continuous functions through the 
thickness. The zigzag kinematic framework enables sufficiently accurate and computationally efficient modeling of 
a wide range of homogeneous and heterogeneous laminates without the use of shear correction factors. Novel zigzag 
functions, derived a priori from constitutive relations without enforcing debilitating stress -equilibrium constraints, 
are responsible for overcoming several critical shortcomings of the earlier zigzag theories. The stress resultants 
obtained from the equilibrium equations are physically consistent with their definitions based on Hooke’s relations. 
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In addition to the in-plane deformations, this element includes both transverse shear and transverse normal 
(thickness-stretch) deformations, and computes accurately the in-plane and transverse stress components through the 
thickness, i.e., it permits the determination of all six stress and strain components. This C°-continuous element 
maintains a fixed number of kinematic unknowns regardless of the number of material layers. Using the principle of 
virtual work, equilibrium equations and consistent boundary conditions are derived in a variationally consistent 
manner. 


As shown in Figure 2, the finite element is a six-node triangle consisting of three corner nodes and three mid- 
side nodes along the edges. At each corner node, eleven degrees of freedom (DOF) are specified. These are the in- 
plane displacements, u u) and v 01 , the transverse displacements, w 10 , vr, (,, ,and w ( f , the bending rotations, 8 , (0 and 
8f , and the zigzag amplitudes (rotations), y/^ , y/^ , y/^ , and y/^ , with i = 1,2,3 denoting the corner nodes. At 
the mid-side nodes, only the transverse displacements, w <,+3) , Wj l+3) , and vv!,' +3) ( i = 1,2, 3 ), are defined. 



Fig. 2. Description of 6-node anisoparametric element with nodal DOF assignment 


The non-homogeneous nature of a laminated composite is considered by including discrete fiber- and resin-rich 
layers. The element is capable of a highly detailed modeling of each ply (fiber- and resin-rich layers regardless of 
the number of plies) as well as the stiffness degradation of any of the layers. Unlike the existing finite elements for 
composite laminates, the present provides robust and accurate prediction of all six stress components (in-plane and 
transverse normal and shear stresses). The key features of this element are summarized as 

• Employs refined zigzag theory 

• Accounts for the through-thickness stretch and shear deformations 

• Uses anisoparametric interpolations that do not suffer from shear locking 

• Requires no shear correction factors to yield accurate predictions 

• Uses zigzag functions that accurately model the layerwise shear stiffnesses 

• Models accurately the fiber- and resin-rich regions, in addition to the highly stiff and very compliant layers 

• Provides highly effective modeling for thick-section laminates, sandwich structures, and bonded joints 

• Models accurately the in-plane, transverse normal, and transverse shear stresses 

• Computes highly accurate interlaminar stresses via integration of equilibrium equations 

Progressive failure criteria 

Failure modes in laminated composite panels are strongly dependent on ply orientation, loading direction and 
panel geometry. There are four basic modes of failure that occur in a laminated composite structure. These failure 
modes are; matrix cracking, fiber-matrix shear failure, fiber failure, and delamination. To simulate damage growth 
accurately, the failure analysis must be able to predict the failure mode in each ply and apply the corresponding 
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reduction in material stiffness. The failure criteria included in the present analyses are those proposed by Hashin [4] 
as 


Matrix compressive failure — (cr 99 + <t 33 ) < 0 

j 2 = I' K f a 22 + P~33 1 , (°~22 + °~33 ) , (°~23 ~ <7 22 (7 33 ) 
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Matrix tensile failure — (<x 22 + <r 33 ) > 0 
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Fiber-matrix-shear compressive failure — cr n <0 
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in which cr with (i,j = 1,2,3) are the stress components in reference to the material coordinate system. The 
parameters (X t , Xf) and (K, Y) represent the tension and compression strengths along the fiber and transverse 
directions, respectively. The in-plane shear and transverse shear strengths are S 12 , S l3 and S 23 , respectively. 

Recursive degradation 

Degradation of ply stiffness is a common approach as part of the laminate failure analysis. Based on the 
stresses and strains in each ply, a lamina failure criterion such as that of Hashin [4] is employed to determine the 
particular ply which fails first and its failure mode. Typical values for the degradation factor can range from a very 
small value ( e.g 10-6) to a large value (e.g., 0.8) for this type of material degradation. As suggested by Knight 
[23], recursive degradation successively degrades the material stiffness coefficients in a gradual manner to avoid 
numerical difficulties arising from abrupt changes in material stiffness. When the degradation factor is not small 
(for example, 0.5), then some of the numerical convergence issues associated with an instantaneous local change in 
material stiffness are avoided. Specifying recursive degradation with a near -zero degradation factor is nearly 
equivalent to specifying instantaneous degradation with the same near -zero factor. 
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Non-linear incremental finite element analysis for stress redistribution 

The laminate with reduced stiffness is again analyzed for stresses and strains. The lamina failure criterion 
predicts the next ply failure, and laminate stiffness is accordingly reduced again. This cycle continues until ultimate 
laminate failure is reached using an incremental approach in accordance with the incremental equilibrium equation 
in the form 

\ Au ( ' +A *' i+1) — JJtf+AO ^ j 


in which K(C“ +a '''’) represents the tangent stiffness of RZE at load step (f + Af) in the i' h iteration with 
Qit+At.o _ qu+m,i) ^E xl ,E 22 ,Ey } ,G l2 ,G i3 ,G 23 ,v i2 ,v n ,v 23 ) being the material stiffness matrix of the k ,h ply in laminate. 

The known external force is represented by R (r+A,) , and the internal force vector by f ,+ ’ . The unknown 
incremental displacement vector at the (7 + 1)* is Au (,+ a, ’ ,+ ’ . Until first ply failure, the material stiffness for each 
ply remains as original and not degraded, i.e., C 2 +v ' n = C[°' . At the onset of ply failure, the sets of ply stiffness 
properties, (E n , E 22 ,G l2 ,G 23 ,G l3 ,v n ) , (E 22 ,v 12 ), (G 12 ,G 13 ,v 12 ) , and ( E 33 j for fiber failure, matrix failure, fiber- 
matrix shear, delamination, respectively, are degraded by a factor of p whose initial value is 0.5. After each load 
increment, these sets of ply properties are recursively degraded as follows 


E% ) =P%E™ with 

with /f =/^ 
vf with 


(5a) 

(5b) 

(5c) 


in which £( 0) , G( 0) and u™ 1 represent the elastic moduli, shear moduli and Poisson’s ratios, respectively, for an 

undamaged ply. The ultimate laminate failure is established when the determinant of tangent stiffness matrix 
becomes nearly zero. 


III. Numerical Results 

The numerical results demonstrate the capability of this approach by first comparing the present predictions 
with the laminate tensile experiments by Sun et al. [18] and the predictions by Triplett [19 ]. Previous experimental 
results by Suemasu et al [20] serve as a validation for a laminate with an open hole under compression. The 
experimental results by Chang and Chang [8] and previous predictions by Sleight [21], and Gunel and Kayral [22] 
are considered in the validation of predictions for a laminate with a hole under tension. 

Strength of laminates under off-axis loading 

The laminates that are considered are made of AS4/ 3501-6 Gr/Epoxy unidirectional plies; they have the ply- 
level properties summarized in Table 1. 


Table 1. Ply-level properties for AS4/3501-6 Gr/Ep system Sun et al. [18] 


E n 

153.7 GPa 


2171.0 MPa 

E 22 

ll.OGPa 

X 

c 

-2013.0 MPa 

G]2 

6.9 GPa 

K 

67.0 MPa 

G2 

0.32 

K 

-206.8 MPa 

ply 

0.13 mm 

s a 

110.3 MPa 


As pointed out by Sun et al. [18], off-axis loading is achieved by rotating all plies of the laminate by the same 
amount as the angle of loading. The finite element model of each laminate is composed of two triangular elements 
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(RZE), and the boundary and loading conditions are shown in Fig. 3. This configuration results in a uniform strain 
field. The predictions from the proposed approach are compared with the actual strength measured by Sun et al. 
[18] and other predictions by Triplett [19], in which only one quadrilateral element subjected to constant strain 
conditions is employed to determine both the first-ply and ultimate failure loads of a laminate; these are summarized 
in Tables 2-5 for various off-axis loading angles. As discussed by Sun et al. [18], the degree of agreement between 
the predicted and measured strength values depends on the off-axis angle. The agreement is remarkably good for 
the majority of the cases, and remains quite respectable for the remaining cases. 



Fig. 3. FE model of a square laminate under uniaxial tension to achieve uniform strain. 


Table 2. Comparison of measured and predicted strength for a laminate lay-up of [0/45/-45/90] 


j Laminate lay-up — [0/45/-45/90] s | 

Off-axis angle 

Ultimate strength, MPa j 

Experiment [18] 

Triplett [19] 

Present 

0 

765 

721 

725 

7.5 

752 

731 

749 

15 

774 

788 

833 

22.5 

832 

885 

987 


S 


Table 3. Co mparison of measured and predicted strength for a laminate lay-up o f [90/30/-30] s 


! Laminate lay-up — [90/30/-30] s | 

Off-axis angle 

Ultimate strength, MPa j 

Experiment [18] 

Triplett [19] 

Present 

0 

966 

1064 

1086 

7.5 

908 

885 

921 

15 

837 

782 

823 

22.5 

807 

731 

776 


Table 4. Co mparison of measured and predicted strength for a laminate lay-up o f [0/45/-30] s 


| Laminate lay-up - [0/45/-30] s | 

Off-axis angle 

Ultimate strength, MPa | 

Experiment [18] 

Triplett [19] 

Present 

0 

883 

756 

756 

7.5 

843 

769 

780 

15 

929 

859 

864 

22.5 

1028 

1026 

1093 

26 

1129 

1167 

1138 

30 

1074 

1026 

1026 

45 

818 

744 

804 
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Table 5. Co mparison of measured and predicted strength for a laminate lay-up o f [90/0/90/0] s 


! Laminate lay-up - [90/0/90/0] s ! 

Off-axis angle 

Ultimate strength, MPa j 

Experiment [18] 

Triplett [19] 

Present 

0 

1126 

1073 

1086 

1.5 

1140 

1038 

1087 

3 

1074 

808 

1089 

4.5 

1018 

600 

1093 

6 

861 

473 

1099 

7.5 

713 

452 

852 

15 

394 

346 

441 

22.5 

288 

288 

312 


Strength of laminates with a hole under compression 

In this demonstration example a composite laminate with a hole under compression is considered and the 
present strength predictions are compared to the experimental strength measurements by Suemasu et al. [20]. The 
laminate is made of T800H/3633 Gr/Epoxy unidirectional plies; the ply-level properties are summarized in Table 6. 
The geometry and finite element model of the specimen are shown in Fig. 4. The width, length, and thickness of the 
laminate are denoted by W = 38 mm , L = 1 1 8 mm , and t = 2.2 mm , respectively. The diameter of the hole located 
at the center is <i =6.35 mm. The laminate lay-up is of [-45/90/45/0| 2s , and it is made of T800H/3633 
Gr/Epoxy. 

As shown in Fig. 4, the specimen is clamped along the left edge and subjected to uniform end shortening of A 
via a rigid clamp along the right edge. Note that the rigid clamp only allows end shortening (horizontal 
displacement) along the right edge. The corresponding load, P , is then calculated as a reaction force in the FE 
analysis. The first ply failure and ultimate strength predictions from the present approach are compared with the 
experimental measurements by Suemasu et al. [20]. In Fig. 5, the present analysis prediction for the quasi-isotropic 
laminate layup is shown with dark solid line. In this figure, the circles represent the load at which the first-ply 
failure occurs; the squares denote the load at which deviation from the linear variation starts; and the triangles 
represent the ultimate load. It is observed that the predicted first ply failure is 16.2 kN and the ultimate strength is 
21.8 kN as compared to 17 kN and 23.2 kN measured by Suemasu et al. [20], respectively. The agreement with the 
experiment is remarkably close. 

In order to further demonstrate the effect of lamination on the strength, two other lay-ups of, [ 90 / 0 / 90 / 0 ] 2l 
and [90/ 30/ -30/ Of, are also considered. The prediction for these layups are shown in Fig. 5 with blue and red 
dashed lines. As expected, these laminates have higher strength than the one tested by Suemasu et al. [20] as shown 
in Fig. 5. 


Table 6. Ply-level properties for T800H/3633 Gr/Epoxy Suemasu et al. [20] 


E n 

148.0 GPa 


2000.0 MPa 

E 22 

9.56 GPa 


-1500.0 MPa 

P 

II 

4.55 GPa 

Y 

50.0 MPa 

g 23 

3.17 GPa 


V \2 = G3 

0.3 

Y 

-150.0 MPa 

^23 

0.49 

I c 

^ ply 

0.1308 mm 

S 12 

100.0 MPa 
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Fig. 4. Geometry, loading, boundary conditions, and finite element discretization for a compressively loaded 

laminate with a hole. 



Fig. 5. Load vs. end-shortening response of compressively loaded laminates with a central hole 


Strength of laminates with a hole under tension 

In this example, a composite laminate with a hole under tension is considered, and the present strength 
predictions are compared with the experimental strength measurements by Chang and Chang [8] as well as previous 
predictions by Sleight [21], and Gunel and Kayral [22], The geometry and finite element model of the specimen are 
shown in Fig. 6. The width, length, and thickness of the laminate are denoted by W = 25.4 mm. , L = 203.2 mm. , 
and t = 2.616 mm. , respectively. The diameter of the hole located at the center is d = 6.35 mm. . The laminate lay- 
up is made of T300/1034-C Gr/Epoxy unidirectional plies with a lay-up of [0/(±45) 3 /90 3 ] s . The ply-level 
properties are given in Table 7. 

While the side (horizontal) edges of the specimen are free from constraints, the loaded edge is subjected to 
uniform end-displacement of A via rigid clamps, and the left vertical edge is clamped as depicted in Fig. 6. As 
shown in Fig. 7, the end shortening-load curve obtained from the present analysis for the laminate layup considered 
previously by others is represented with dark solid line. Similar to the previous example, the first ply, first 
deviation, and ultimate loads are, respectively, denoted by circles, squares, and triangles in Fig. 7. 
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The first ply failure and ultimate load predictions from the present approach are compared with the 
experimental measurements by Chang and Chang [8] and previous predictions by Sleight [21] and Gunel and Kayral 
[22], As shown in Fig. 7, the predicted first ply failure is 10.72 kN and the ultimate strength is 16.7 kN as compared 
to the measured ultimate strength of 15.6 kN (experimental first ply failure is not available). The previous 
predictions by Sleight [21] and Gunel and Kayral [22] are, respectively, given as (6.7kN, 14.3kN) and (8.0kN, 
16.5kN), in which the first number in the parenthesis denotes the first ply failure load and the second number 
represents the ultimate load. It is observed that the comparisons against the experiment and other predictions 
indicate reasonably close agreement. 

In order to further demonstrate the effect of lamination on the strength, two other lay-ups of, [0, / (±45) , /90, [ 
and [0, /(±45), / 90] y are considered. In Fig. 7, the results for these two laminate layups are represented by dashed 
red and blue lines, respectively. As expected, these laminates have higher stiffness than the one tested by Chang and 
Chang [8], 


Table 7. Ply-level 

properties for of T300/1034-C Gr/Epoxy, Chang and Chang [8] 

E n 

146.85 GPa 

X, 

1723.7 MPa 

E 22 

1 1 .4 GPa 


-1378.9 MPa 

r Q) 

II 

6.2 GPa 

v 

66.5 MPa 

G 23 

2.18 GPa 


G2 = Pl3 

0.3 

Y 

-268.2 MPa 

^23 

0.49 


^ ply 

0.1308 mm 

S n 

133.7 MPa 



Fig. 6. Geometry, loading, boundary conditions, and finite element discretization for a laminate with a hole 

under tension. 
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Fig. 7 Load vs. end shortening response of laminates with a central hole under tension 


IV. Concluding remarks 

The refined zigzag element, RZE 12 21 , accounts for the discrete nature of fiber- and resin-rich layers of each ply 
as well as the stiffness degradation of any of the layers. This new element predicts accurate in-plane, transverse- 
normal, and transverse-shear stresses; thus, it eliminates the lack of reliability associated with inaccurate 
computation of stress redistribution while performing progressive failure in composite laminates. This approach will 
enable the designers to consider various laminate lay-ups without relying on significant amount of material testing 
and characterization. As a result, significant reduction in cost and design cycle time are expected in the use of high- 
performance material systems as they become commercially available. 
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